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Abstract 

The asymptotic behavior of open plane sections of triply periodic sur- 
faces is dictated, for an open dense set of plane directions, by an integer 
second homology class of the three-torus. The dependence of this ho- 
mology class on the direction can have a rather rich structure, leading in 
special cases to a fractal. In this paper we present in detail the results 
for the skew polyhedron {4, 6 | 4} and in particular we show that in this 
case a fractal arises and that such a fractal can be generated through an 
elementary algorithm, which in turn allows us to verify for this case a 
conjecture of S.P.Novikov that such fractals have zero measure. 

1 Introduction 

The study of plane sections of triply periodic surfaces in R 3 was initiated by 
S.P.Novikov in Nov82] who raised a question whether open (unbounded) compo- 
nents of such a section have some nice asymptotic behavior. This was motivated 
by an application to conductivity theory. A number of general theoretical results 
has been obtained since then by A.Zorich |Zor84] . I. Dynnikov |Dyn97| [Dyn99| , 
and R.DeLeo DeL04, DeL05]. For a number of surfaces, R.DeLeo performed 
a numerical simulation, which confirmed the general conclusions of the the- 
ory [DeL03j. 

In |DeL06j the main results were generalized to polyhedra. Among the class 
of piecewise linear triply periodic closed surfaces, the one of infinite skew poly- 
hedra |Cox37j is the most suitable for numerical explorations of the geometry 
of plane sections. In [DeL06j the case of the regular skew polyhedron V of type 
{6, 4 | 4} was studied numerically in detail, showing that the dependence of open 
section's asymptotics on the plane direction keeps its rich fractal-like structure 
also in the piecewise linear case. 
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In this work we present the results for the dual of V, namely the cubic 
polyhedron C = {4, 6 | 4} |CGS03j ; note that both V and C are rough PL- 
approximations of the smooth surface {cos a; 1 +cosa; 2 +cosa; 3 = 0}, which was 
itself studied numerically in |DeL03j . It turns out that the C case is rather note- 
worthy because its correspondent fractal can be generated recursively through 
a simple algorithm, which, on one hand, allowed us, for the first time, to ver- 
ify, in this concrete case, the conjecture [NM03] that the set of exceptional 
directions has zero measure, and, on the other hand, made possible a system- 
atic comparison with the numerical data obtained through the NTC software 
library |DeL99j . 

2 Topological structure of plane sections of triply 
periodic surfaces 

Let M C T 3 be an embedded closed null homologous surface in the torus T 3 , 
H = [hi, /12, ha) G (M 3 )* a covector. We denote by M the preimage of M 
under the projection R 3 — > T 3 = R 3 /Z 3 . We consider the sections of M by 
planes (H,x) = const (we call them H -sections) and we are interested in the 
asymptotical behavior of their unbounded regular connected components (if 
any). Since only the direction of covector H matters, sometimes we shall treat 
H as a point of the projective plane RP 2 . 

In studying this question, the foliation Tr induced on M by the closed one- 
form ui = (hidx 1 + h 2 dx 2 + h 3 dx 3 )\M plays the crucial role. It is well known 
that, with probability 1, in a proper sense, a smooth closed one-form whose 
critical points are all saddles induce dense leaves on M. However, by restricting 
attention to a special class of one-forms that are pull-backs of a constant one- 
form on T 3 , we fall exactly in the opposite situation, namely, with probability 
1, open leaves are either absent or confined to genus one minimal components 
of the foliation on M, and dense leaves arise only in exceptional cases. 

There are three principally different types of foliations Tu and corresponding 
-ff-sections that may arise, which we call trivial, integrable, and chaotic. Most 
typically, trivial means that all regular leaves of Fjj are closed, integrable means 
that minimal components of Th filled with open leaves are of genus one, and 
in the chaotic case there is a minimal component of genus > 1. More precise 
definitions are as follows. 

Let N C T 3 be a piece-wise smooth embedded surface in T 3 such that 
N\M consists of disjoint open disks each of which lies in a plane of the form 
{x e M 3 ; (H, x) = const}. Such a surface is obtained by, first, cutting M along 
some closed null homologous leaves of Th or null homologous saddle connection 
cycles, second, removing some of the obtained connected components, and then 
gluing up planar disks in order to obtain a closed surface. For such a surface 
N, any leaf of Th is either contained in N or disjoint from N. In the former 
case we say that the leaf is absorbed by N. When saying this we shall assume 
that N is of the just specified form. 
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Trivial case. Every leaf of Th is absorbed by a two-sphere or a null homolo- 
gous two-torus. If covector H is totally irrational, i.e., dirriQ(/ii, /12, /13) = 
3, then this just means that all connected components of all ii-sections of 
M are compact. If dimQ(/ii, /12, /13) < 3, then, additionally, periodic, i.e., 
invariant under a non-trivial shift, unbounded component of iJ-sections 
may arise; 

Integrable case. Every leaf of Th is absorbed by a sphere or a two-torus, 
and at least one leaf is absorbed by a two-torus with non-zero homology 
class. In this case, every regular non-closed component of an H -section is 
a finitely deformed straight line, i.e., it has the form j(t) = t ■ v + O(l) for 
some parametrization, where v is a non-zero vector. If dimQ(/ii, hi, /13) = 
3, then the non-zero homology class of the tori absorbing leaves of J-r 
is uniquely defined up to sign. We denote it by Lm,h and consider as 
an integral covector in M 3 . The identification of F 2 (T 3 ,R) and (K 3 )* = 
i/ 1 (T 3 ,R) is given by the Poincare duality. This covector must obviously 
vanish at vector v: (Lm.h,v) — 0. So, if we assume our three-space 
Euclidean, then (unless H and Lm,h are colinear) we can simply write 
v = Lm.h x H. If dirriQ(/ii, /i2, /13) < 3 the covector Lm,h may not be 
uniquely defined (up to sign), but there may be at most two different 
choices. We denote the projective class (L\ : L 2 : L 3 ) 6 QP 2 of Lm,h by 
@m,h and call the soul of the foliation !Fh . 

Chaotic case. None of the above. If dimq(hi, /12, /13) = 3, this means that 
a minimal component of Th has genus ^ 3. The behavior of the cor- 
responding _ff-sections have not been studied, but the known examples 
suggest that, typically, a chaotic ii-section contains a single unbounded 
curve that "wonders all around the plane", i.e., a ^-neighborhood of the 
curve is the whole plane for some finite d. 

For a fixed surface M and a rational point I S QP 2 C MP 2 we denote by 
T>M,e the set 

V M .i = {He KP 2 ; £ m ,h = £}. 

If ^m,h is not uniquely defined then the point (hi : li2 : ft.3) is attributed to both 
corresponding subsets. The set of points (hi : h% : h%) such that the -ff-sections 
of M are chaotic will be denoted by £(M). 

The following three propositions are extracted from |Dyn99| . 

Proposition 1. For a generic surface M C T 3 the sets T>m,i cLre disjoint closed 
domains with piece-wise smooth boundary. The set £(M) is disjoint from QP 2 
and has zero measure. The set of directions H with trivial H -sections is open. 

In other words, the first claim of this proposition says that t, where defined, is 
a locally constant function of H. We call the non-empty domains "Dm,£ stability 
zones and refer to £ as the label of the stability zone T>m,i- 

For studying the stability zones, it is usefull to consider a 1-parametric family 
M c = {x 6 T 3 ; f(x) = c] of level surfaces of a fixed smooth function. 
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Proposition 2. For a generic function f, there are continuous functions e\,ei ■ 
RP 2 -> R such that 

• e x {H) < e 2 (H) for all H e RP 2 ; 

• the H -sections of M c are trivial if and only if c ^ [e\{H) 1 e^iH)]; 

• if ei(H) < e2(H), then the H-sections of M c are integrable for all c € 
[e\{H) 1 62 (i?)], and the soul £ of the corresponding foliation T C ,H is inde- 
pendent of c. 

Wc define generalized stability zones T>f/ as T>f t i — U c 2?m c /, and the set 
£{f) as £{f) = U C £(M C ). 

Proposition 3. For a generic f , generalized stability zones are closed domains 
with piece-wise smooth boundary. If £ ^ £' then the zones T>fj T^ji' can 
only have intersections at the boundary, and, moreover, the number of their 
common points is at most countable. If the whole RP 2 is not covered by a single 
generalized stability zone, then the number of zones is countably infinite, and 
the set £(f) = MP 2 \ (UeT>f t i) is non-empty and uncountable. 

It may happen that there is just one generalized stability zone (say, a small 
enough perturbation of the function sin(a; 1 ) will work), but it is also easy to 
find a function / with non-empty <?(/): any function with cubical symmetry is 
such. In all examples known to us two different generalized stability zones have 
at most one point in common. 

It follows from Propositions [T] [3] that the union U£ g Qp2int(I?/ i £) of the interi- 
ors of the zones is an open everywhere dense subset of RP 2 and its complement 
£{f) has the form of a two-dimensional cut out fractal set. 

Pro posit ion 4 (DcL04j). If th ere is more than one generalized stability zone, 
then £{f) is the set of accumulation points of the set of their souls. 

It is plausible but still unknown whether £(f) has always zero measure. The 
following stronger conjecture was proposed in [NM03] . 

Conjecture. Whenever £(f) is non-empty, the Hausdorff dimension of £(f) 
is strictly between 1 and 2 for every f. 

So far only numerical checks of this conjecture were available in the litera- 
ture |DeL03[ [DcL06j: in the next sections we shall provide, for the particular 
case of the polyhedron C, a full proof of the weaker, zero measure, conjecture 
and good numerical evidence for the stronger one. 

3 Stability zones of C 

The regular skew polyhedron C = {4, 6 | 4} (see Figure[lJ is, up to isometries, the 
unique cubic polyhedron with all monkey-saddle vertices |CGS03| : the vertices 
of its embedding in T 3 = [0, l] 3 / ~ are the eight points in the orbit of P = 
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Figure 1: The {4, 6 | 4} polyhedron embedded in the three-torus. Shown is the fun- 
damental domain, which lies in the [0, l] 3 cube. 

(1/4,1/4,1/4) under the cubic symmetry group. As level surface, C can be 
represented in the [0, l] 3 cube as _1 (O) for 

9(x 1 7 x 2 , x 3 ) = mid(|2x 1 - 1|, \2x 2 - 1|, \2x 3 - 1|) - -, 

where mid(a, b, c) is the middle value among a, b, and c. 

This surface, as well as the surface cos x 1 + cos x 2 + cos x 3 — and C 's dual — 
the truncated octahedron, has a very strong symmetry, namely, its exterior is 
equal to its interior modulo a translation. This means that for 9 the functions 
e± t 2 mentioned in Propositoin [2] are such that e% = — e%, and hence, stability 
zones of the surface C coincide with generalized stability zones of the function 
0. 

Let us denote by tpi , "02 , "03 the following projective transformations: 

ipi(hi : h 2 : h 3 ) = (hi : h 2 + hi : h 3 + hi), 
ip 2 {hi : h 2 : h 3 ) = (hi +h 2 :h 2 :h 3 + h 2 ), 
ip 3 (hi : h 2 : h 3 ) = (hi + h 3 : h 2 + h 3 : h 3 ). 

Theorem 1. For the surface C the stability zones are as follows: 

£>(i:0:0)(C) = {(hi : h 2 : h 3 ) £ RP 2 ; hi > \h 2 \ + \h 3 \}, 

(C) = {{hi : h 2 : ha) € RP 2 ; \hi\ + \h 2 \ + \h 3 \ < ihi, 4h 2 , 4/i 3 }, 

©Vn(V-n(-^ fc (( 1:1:1 ))-))( C ) = V'ti (^tiG'-V't* (©(1:1:1) (C)) •••)): 
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Figure 2: A connected component of C after cutting along closed leaves 

where . . . , i^) is an arbitrary finite sequence of elements from {1, 2, 3}, and, 
in addition, all zones obtained from the listed ones by cubical symmetries: per- 
mutations and changing signs of the coordinates. 

The proof of this theorem will rest on the following two lemmas. 
Lemma 1. We have 

Z>(i:0:0)(C) 3 {(/ii : h 2 : h 3 ) e MP 2 ; h x > \h 2 \ + \h 3 \}. 

Proof. If the inequality hi ^ \h 2 \ + \h 3 \ is satisfied then the plane hix 1 + h 2 x 2 + 
h 3 x 3 — const passing through the center of the unit cube [0, l] 3 separates the 
faces x 1 = and x 1 = 1. This means that closed leaves of the corresponding 
foliation will cut our triply periodic surface C into parts each of which is a finitely 
deformed plane x 1 = const with holes (see Figure [2]). Filling the holes by flat 
disks and projecting to T 3 we obtain two tori whose homology class is equal up 
to sign to (1,0,0) € (M 3 )*. □ 

Lemma 2. Let h\,h 2 ,h 3 ^ 0, (hi,h 2 ,h 3 ) 7^ (0,0,0). Then we have (hi : h 2 : 
h 3 ) £ T>t if and only if ipi((hi : h 2 : /13)) S V^m, where i = 1,2,3. 

Proof. For convenience, we shift the coordinate system as follows: (x , x 2 ,x 3 ) 1— > 
(x 1 + 1/4, a; 2 + 1/4, x 3 + 1/4). Our surface C cuts R 3 into two parts, iV_ and 
N +1 that can now be characterized as follows: iV_ (resp. N + ) consists of points 
(x 1 , x 2 ,x 3 ) £ M 3 such that at least two (resp. at most one) of the three numbers 
{x 1 }, {x 2 }, {x 3 } are in the interval [0,1/2] (resp. (0,1/2)), where {x} denotes 
the fractional part of x. 
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We may assume i — 3, (hi, h 2 , h 3 ) = (a, (3, 1) without loss of generality. Let 
II be the plane defined by ax 1 + (3x 2 + x 3 = c. Denote by Q- the projection of 
LT n iV_ to the x 1 , x 2 -plane along x 3 . According to the description of N_ given 
above Q- is the set of points (x 1 ^ 2 ) £ R 2 such that exactly two of the three 
numbers {x 1 }, {x 2 }, {c — ax 1 — f3x 2 } are in the interval [0, 1/2]. 

Denote by the square {(a; 1 ,^ 2 ) £ M. 2 ; a < x 1 ^ a + 1/2, b ^ x 2 < 
b + 1/2}, and by S m the strip defined by m c — ax 1 — fix 2 ^ m + 1/2. We 
have 

Q- = ( (J Oi,*) u ( (J □ J +i/ 2 ,fc n s m ) u ( (J n.-fe+i/a n sjj. 

The first part in this union, (J . fcgZ does not depend on II. We call these 
squares mainlands. 

Each intersection □j +1 / 2 ,fc H with j, fc, m € Z, whenever non-empty, is a 
convex polygon that can be of the following three types: 

cape: it has a piece of boundary in common with exactly one of the mainlands 
or □j+i,k; 

bridge: it has a piece of boundary in common with both mainlands LTj- ^ and 
island: it is disjoint from mainlands. 

Similarly, one defines the type of a polygon Oj^+1/2 H 5 m regarding adjacency 
to the mainlands 0,-^ and see Figured] 

Capes are not interesting for us because their removal is equivalent to a finite 
deformation of Q_. It is not hard to write the necessary and sufficient condition 
for □ :)+1 /2,fc H S m to be a bridge: 

c - a(j + ~) - p(k + ~) - i < to < c - a(j + 1) - /3fc, (1) 
or an island: 

c - «(i + 5) - /3(fc + ~) - 5 > to > c - a(j + 1) - (2) 

Now we apply i/'s to H, which gives H' = (a', /?', 1) = (a + 1, /3 + 1, 1). Let 
Q'_ and S' m be defined in the same way as Q- an d S m with a and /? replaced 
by a' = a + 1 and (3' = [3 + 1, respectively. Then the intersection Oj +1 / 2 ,k H 
is a bridge if and only if 

c-(a+l)(j + i)-C9+l)(fc+i)-i <m<c-(a + l)(i + l)-(/3+l)fc, 
which can be rewritten as 

c-a(j + -) - /?(fc + - - m + j + k + 1 sC c- a(j + 1) - f3k. 
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i i i i 

i i i i 

i i i i 

3 3+1 3+1 3+1 

Figure 3: Geography of the section II n AT_ 

Thus, □j+i/2,fenS' m is a bridge if and only if so is ^j+i/2,k^S' m -j-k-v Similarly, 
the same is true about islands as well as bridges and islands in squares of the 
form Of^+i/a. 

So, bridges and islands of Q'_ in the square □ J +i/2.fc or are in a 

natural one to one correspondence with those of Q-. Therefore, Q'_ and Q_ are 
obtained from each other by a finite deformation. The geometrical difference 
between Q'_ and Q_ can be vaguely described as follows: islands and bridges 
of Q'_ are narrower than those of Q-, and Q'_ has more capes. See FigureHJfor 
an example. 

In the genus three case the integrability of our foliation is equivalent to the 
existence of closed fibres of the foliation (or null homologous saddle connection 
cycles). We have just seen that 77-sections and H '-sections are obtained from 
each other by a finite deformation. Hence, they are both integrable or both 
chaotic. If they are integrable, let I and I' be the labels of the corresponding 
zones. We want to show that £' = ip^,(t). Since both t and £' are locally 
constant functions of H it is enough to consider the case of totally irrational 
H. Then the asymptotic direction v = H x t is of irrationality degree two, i.e., 
diniQ^ 1 , v 2 , v 3 } = 2, and I is the only rational covector up to multiple that 
vanishes at v. Thus, it suffices to show that ^3^) vanishes at v', the direction 
of open components of iJ'-sections. 

The latter follows easily from the fact that the projections of v and v 1 to the 
x , x 2 -plane coincide, which, in turn, follows from the coincidence, up to finite 
deformation, of Q_ and Q'_. □ 
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Q- 



Figure 4: The transition from Q_ to Q'_ produces more capes and makes the 
bridges narrower 

Proof of Theorem [7J Due to the cubical symmetry of the surface it suffices to es- 
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1)} and i is not of the 
form . .ip ik ((l : 1 : 1))...)). We want to show that !Z e = Vg for 

all £ £ QP 2 n C + . We have already shown in Lemma Q] that TZg C T>g for 
£ = (1 : : 0) By symmetry it is also true for £ = (0 : 1 : 0), (0 : : 1). 
It is a straightforward check to see that 

ft(lilil) = ^1 (^(1:0:0) n C+) U ^(^(0:1:0) H C+) U ^(^(0:0:1) H C+). 

By Lemma[5]this implies 7^(i : i : i) C X > (i : i : i) as (1 : 1 : 1) = V>i(l : : 0) = ^(O : 
1:0) = ^3(0 : : 1), and, therefore, K t C V e for all feC+fl QP 2 . 

In order to establish Theorem Q] it suffices to show that the zones T>t are not 
larger that 7Zg, and there are no other stability zones. Both claims follow from 
the fact that [J e TZe covers all rational points: 



\jnt d c+n< 



Indeed, let cp be the following map from C+ to itself: 

(hi : hi — hi : h 3 — hi), if < hi < h 2 , h 3 , 
(h 2 : h 3 : h\), otherwise. 



<p(hi : h 2 : h 3 ) 
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Figure 5: Picture of the fractal in the disc model of RP 2 . The center of the disc 
corresponds to the z axis, so the central square is the stability zone Z>(o : oa)j t ne one 
touching it in the left and right vertices is 23(i : o : o) and the third one is 2?(o : i : o)- The 
green triangles are the stability zones 2?(i : i : i), X>(i:i : _i), 25(i : -i : i) and 25(i : -i : -i). 



By construction, for every H e C+ we have H G if and only if ip(H) g 

U e TZi- If H is a rational covector from C+, then after applying ip finitely many 
times, one of the coordinates of the obtained covector becomes zero. All such 
points are covered by T^^o), 7^(0:1:0) and 7?.(0:0:i)- D 

So, we have the following picture in MP 2 : four lines hi ± /12 ± /13 = cut 
RP 2 into three "squares", which are zones 2?(i : o : o)7 ^(o-.i-.o), ^-'(0:0:1)7 an d four 
triangles obtained from each other by cubical symmetries, in which there are 
infinitely many stability zones. We shall concentrate on the triangle that is 
contained in C+. This triangle has vertices (1 : 1 : 0), (1 : : 1), (0 : 1 : 1) and 
is defined by the inequalities 

— — I ->h u h 2 ,h 3 >0. 

Let us denote this triangle by A. The zone is also a triangle that is 

contained in A and has its vertices, (2 : 1 : 1), (1 : 2 : 1), (1 : 1 : 2), at the sides 
of A. The complement A \ Pn : i : i) consists of three triangles that are exactly 
^i(A), V 2 (A), and </> 3 (A). In each triangle Ai = Vi(A), A 2 = ?/> 2 (A), A 3 = 
^3 (A) the picture is obtained from that in A by the corresponding projective 
transformation ipi. 

For a finite sequence a = (oi, a>2i ■ ■ ■ > a>k) of indices 1,2,3 we denote by 
ip a the mapping ip ai o ip a2 o . . . o ip ah . By a', a", and a'" we denote the se- 
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quences (aj., . . . , ajt, 1), (oi, . . . ,a&, 2), (oi, . . . , cifc, 3), respectively. For any such 
sequence a we have the following. The triangle A a = ip a (A) is bounded by 
the zones Dtp a (i:0:0); 2?Va (0:1:0); 2?Va (0:0:1) • It contains the zone 2?^> (i:i:i) whose 
vertices ^ a (2 : 1 : 1), ^> (1 : 2 : 1), "0a(l : 1 : 2) are on the sides of A a , and the 
complement A a \ £\/, a (i : i : i) is the union of the triangles A /, A a », A a >/>. 

Proposition 5. The intersection £{C) n consists of points of the form 
lim ^(^(...^((1 : 1 : 1))...)), 

A: — >oo 

where (ii,«2; ■ ■ ■) fwis over a// possible sequences of elements from {1, 2, 3} con- 
taining each index infinitely many times. Other points in £(C) are obtained from 
these by cubical simmetries. 

Proof. From the structure of stability zones established above it follows that 
the intersection £ n C+ is the union of subsets 

f|V«(V»fa(...^(A)...)) 
k 

over all sequences (i^) G 3 N in which all three indices appear infinitely many 
times. It suffices to show that every such a subset is actually a single point, 
which follows from Proposition [5] below. □ 
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3.1 Measure of S 

In all cases studied so far no algorithm was found to generate all stability zones 
and nothing could be said about the measure of the £(/). This is therefore the 
first case in which it is possible to check the truth of the zero measure conjecture. 

Theorem 2. The set £(C) of "chaotic" directions has measure zero. 

Proof. Again, due to the symmetry, it suffices to prove the claim for E n C+. 
We denote this set by £ + . As we have seen, it is contained in the triangle A, 
and the following holds: 

Applying this recursively, we get 



8+ = M£+) u M£+) u V>si(£+) u ?/> 32 (£+) u V>33i(£+) u ^(s+P 

■ ■■u V^3i(^+)u V^3 2 (f+)u V3^(f+)- (3) 

k fc fc+1 

Let // be the measure of £ + and the measure of V3...3 which is 

k 

equal to the measure of V3...3 2{£+)- The measure of the triangle V'3...3(^+) 

k fc+1 

tends to zero as k goes to 00. Therefore, we have 

00 

k=0 

The idea now is to show that n k ^ c k [i, where Co, c\, . . . are constants such that 



1 

J2 ck < 9' 



k=0 

which, together with the previous equality, implies /1 = 0. 

Now we provide the necessary technical details. First of all, we need to 
parametrise the triangle A. We use the following parametrization: (it, v) h+ 
(1 — v : 1 — u : u + v), u,w>0, u + v ^ 1. 

The property of a set to be of zero measure does not depend on the choice 
of a regular Borel measure. In order to get the desired inequalities we use a 
somewhat artificial measure on A, namely the following one: 

du dv 



(1 + u + v) 2 ' 



By doing so we keep the symmetry between ipi and ip2 (one is conjugate to 
the other by the permutation u ^+ v), so the measures of the sets V3...3 i(£+) 
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and V3...3 2(£+) coincide. We denote by fk the mapping V3...3 1 R, where 
R(hi : fi2 : h 3 ) = (h 3 : hi : h 2 ), written in the u, v-parametrization: 

/fc(u ' U) = L+(fc + l)(v + l)' u + (k + l){v + 1)) ' 

Since we have R{£) — £, the image fk{£+) coincides with V3...3 i(£+)- 

fc 

Obviously, the measure of fk(£+) is bounded from above by CkH, where 
c fc = max J fh {u,v)- 

\ l + « + (fe+l)(D + l) + «+(fc+l)(t> + l) J 

By Jf we denote the Jacobian of the mapping /. A routine check gives: 

(l + u + vf 



Ck = max 



v^o,u+v^i (u+(k + 2)(v + 1)) 2 (m +{k + l)(v + 1)) 

1 



( (fc + 2)(fc + 3) 



2 



if k = 0, 
otherwise, 



• 1 



253 2 2 „ . .„ 1 



tt z w 0.448 < 



4 ^ (fc + 2 )(fc + 3)2 36 3 ~ 2' 
which completes the proof. 



□ 



3.2 Asymptotics of stability zones and fractal dimension 
of£ 

From what said at the beginning of the section it follows immediately that all 
stability zones, except the biggest ones, which are squares, are triangles and 
their labels satisfy a simple recursive relation. 

Proposition 6. Every stability zone ofC, except for the three squares with souls 
(1 : : 0), (0 : 1 : 0), and (0:0: 1), are triangles; these triangles are generated 
starting from A ( and its three symmetric triangles with respect to the coordinate 
planes) according to the the following recursive algorithm: 

1. evaluate the vector sums w\ = v 2 + V3, w 2 = V3 + V\, W3 = V\ + v 2 of all 
pairs of the three vertices Vi of A; 

2. consider the Wi as projective coordinates and add the triangle having those 
points as vertices to the list of all stability zones - its soul is given by the vector 
sum of the three vertices of A; 
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0.40 0.45 0.50 0.55 0.60 

Figure 7: Detail, in the z = 1 projective chart, of the first few zones of the Tribonacci 
sequence of stability zones starting by X>(i : o : o), 2?(0:i:0) an d ^(o-.o-.i)- The labels of the 
zones shown above are, in decreasing area order, (1 : 2 : 2), (2 : 3 : 4), (4 : 6 : 7), 
(7:11: 13) and (13 : 20 : 24). All centers of the zones of the sequence lie on the smooth 
"Tribonacci projective spiral" drawn above which is converging to (a 2 — a — 1, a — 1) ~ 
(.543,.839). 



3. consider now the three triangles Ai = (v\, W2, W3), A 2 — (u>i, V2, W3), A3 = 
(wi, W2, V3) and repeat the algorithm for each of them. 

This fact can be exploited to find an explicit expression for elements of £ 
by considering "spiralling down" (or "steepest descent") sequences of stability 
zones (see Figure [7]). Select indeed an ordered triple (a, /3, 7) of stability zones 
which bind a triangle, namely such that two touch each other and the third 
touches both, and build out of them the recursive sequence of stability zones 
such that every new stability zone is the one whose vertices touch the previous 
three stability zones. It is easy to see that the sequence of these stability zones 
are arranged in a sort of spiral and by construction, the label of every stability 
zone of this sequence is the sum of the labels of the previous three zones, so 
that the sequence of the labels is a Tribonacci sequence in P£z2(T 3 ,Z). 

Proposition 7. The limit of every such sequences belongs to E. In particular 
(a 2 — a — 1 : a — 1 : 1) £ f , where a is the Tribonacci constant, namely the real 
solution of the Tribonacci equation x 3 — x 1 — x ~ 1 = 0. 

Proof. All vertices of stability zones of C are 1-rational points of MP 2 and 
therefore points on the boundaries have irrationality degree not higher than 
2. Now, consider the sequence starting from the three squares a\ = £>(i : o:0)5 
<22 = 2^(0:1:0) i a 3 = 2?(0:0:1)j s0 that the first few next terms of the sequence will 
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be d4 = as = 2?(i : 2:2)) a 6 = ^(2:3:4) an d so on: a simple calculation 

show that the label of the nth stability zone of the sequence, modulo terms in 
n and n , will be 

( §1 _ a n. -13-0 _ n . 1 



\a-0){a-0) (a-0)(a-0) (a-0)(a-0) 



where a = (1+ y/l9 - 3\/33+ \/l9 + 3\/33)/3 ~ 1.84 is the Tribonacci constant 
and 0, the remaining two solutions of the Tribonacci equation x 3 — x 2 — x — 1 = 
0. Since \0\ < a the sequence of labels converges to the triple of coefficients of 
a n , namely = (00 : —0 — 0:1) = (a 2 — a — 1 : a — 1 : 1); these three 
coordiantes are rationally independent so that has rationality degree 3 and 
therefore it cannot belong to any boundary and for Proposition[4]it must belong 
to £. □ 

Note that, since £ is invariant by the ip a , this way we automatically get a 
countable number of explicit elements of £ . 

Proposition 8. Let {2?i}i S N be the set of all stability zones sorted according 
to any recursive algorithm, e.g. Da 1 ...i n ) 3 = ^j„+i o • • • o V , -ti+l(2'(i:i:i))) where 
(ii . . . i n )s is the base 3 expression for the index of the stability zone, and be £ n 
the label associated to V n . Then 21og3(l + 2n) + 1 < ||4|| 2 < 3(1 + 2n) 2lo ^ a 
where a is the Tribonacci constant. 

Proof. A simple induction argument shows that, at every recursion level k, the 
biggest label belongs, modulo permutations of the projective coordinates, to 
the £:th stability zone of the steepest descent sequence having as first three 
elements -D(o,o,i); -D(o,i,o) an d -Dn^o) • Since at the fcth recursion level there 
are 3 fc stability zones it follows that ||£ 3 fe_i || ^ V3a k and therefore \\£ n \\ ^ 

The slowest sequence which can be formed by picking a stability zone for 
every recursion level is instead the one where D/. is the spawn of, say, -D( 0) i), 
15(0,1.0) an d 2?k-i. In this case indeed ||£ 3 fc_i || = i/2(fc + l) 2 + 1 from which 

2 

follows immediately the left part of the inequality. □ 

Proposition 9. Be T>i a stability zone with label t, pi its perimeter and ai its 
area in RP 2 . Then there exist four constants A, B, C , D such that 

A B C D 



\i\\ ' IKII 3 IKII 3 

Proof. The inequalities concerning pg are already known to be true in the general 
case. 

To prove the ones relative to the area we change coordinates in RP 2 so that 
the three points (1 : : 1), (0 : 1 : 1), (1 : 1 : 0) become (0 : : 1), (1 : : 1), 
(0:1:1). This way we can use the Lebesgue measure instead of the projective 
one and we can use the fact that ^5 x, y ^ z; now, be (xi : yi : zA, i = 0, 1, 2, 
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10 20 30 40 50 2 4 6 8 10 12 



Figure 8: (a) Plot of d n = 2 — log 12 V(£ 1 . 2 -n) against n, where V{£ 12 -n) is the 
volume of the neighborhood of £ of radius 1.2~"; for "well-behaving" fractals this 
sequence converges to their Minkowsky dimension. The horizontal line shown in the 
picture has height 1.72. (b) Plot of log 2 N„ against n, where N n is the number of 
squares of side 2 _n which are not completely inside one of the 797161 stability zones 
found at the 12th recursion level. The angular coefficient of the straight line fitting 
the data in the picture is again 1.72. 



the projective irreducible coordinates of the vertices Ai of the triangle T which 
encloses T>£, so that I = Q^j^i : YliVi ■ Si the areas of T and T>i are 
then easily computed respectively as A(T) = Ilf =0 j- and A(T>e) = Ilf =0 - . + * 
(where the indices are meant modulo 3). 

It is convenient for our purposes to consider the maximum norm \\£\\oo since 
in the region under consideration ^ x,y ^ z, so that we can assume without 
loss of generality that Xi, yt ^ Zi, i = 0, 1, 2, and therefore we get trivially that 

) - 1 > 1 - 1 

(z + zi)(zi + z 2 )(z 2 + Zo) ~ (zq + Zl + z 2 ) 3 \\£\\lo 

The second part of the inequality comes from the fact that if max Zi ^ mid Zi + 
min Zi holds for the first triangle of the algoritm in Prop. [6] then it holds for 
every other triangle generated by the algorithm. Indeed assume to fix the ideas 
that zq ^ zi ^ z 2 ] then the new z coordinates of the three vertices under the 
transform ipo( z o '■ z i '■ z i) — (^o : z i + z o ■ z 2 + zq) satisfiy trivially again the 
relation z' 2 ^ z' + z[, and the very same happens for the tranforms under ipi. 
In case of ip2 we have {z' — zq + z 2 , z[ = z\ + z 2 , z' 2 = z 2 } so that the new larger 
z is now z[ but again z[ ^ z' 2 + z' a . Then finally 

M - =(l + (l + -^-) (l + < 8 



(Zq + Zl)(zi + Z 2 )(z 2 + Zq) V Z 1 +Z 2 J\ Z 2 + Z Q J\ Zq + Z 1 

so that A(T>n) ^ 8/pH^. The final statement now follows from the fact that 
all norms are equivalent in finite dimension. □ 

We could not find any way to evaluate exact non-trivial bounds for the 
fractal dimension of £ but numerical calculations suggest that the dimension be 
smaller than 1.8. 
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A simple way to evaluate numerically fratal dimensions is using the Minkowsky 
dimension, namely the limit 

diniM £ — 2 — lim — - — — — 

e-^o log e 

where V(£ e ) is the surface of the e neighborhood of £. In order to do that we 
use the fact that, if A is the area of A and p its perimeter, then 

V(£ e )=pe + eJ2p 1 +A-J2 a * + t 2 ^-J2§:) 

i—1 i—1 i—1 

where k e is the integer such that Pk e +i ^ e ^ Pk e and pk is the radius of 
the inscribed circle to the triangle T>k- In order to avoid infinities we make a 
projective change of coordinates so that the triangle A has vertices in (0 : : 1), 
(1:0:1) and (0:1:1). In Figure [5£a) we show the numerical results we 
got by evaluating the volume of the neighborhoods of £ of radii r n — 1.2 - ™ for 
n = 1, . . . , 50, which suggests a fractal dimension between 1.7 and 1.8. 

A second simple way is to evaluate an upper bound for the box- counting 
dimension, namely the limit 

dim b £ = lim — - — — — - 

e^o — log e 

where N e (£) is the smallest number of squares of side e needed to cover £ (even 
in this case we make the same change of coordinates to avoid infinities). In 
Figure EJb) we show the results relative to covering with squares of side 2~™, 
n = 1, • ■ • ,12, the complement of the 797161 triangles obtained by applying 12 
times the recursion algorithm. Again we get a clear indication of the fractal 
dimension to be between 1.7 and 1.8. 

4 Numerical generation of stability zones 

As an important byproduct of the study of stability zones of C, we were able for 
the first time to compare very accurately the results of the NTC code |DeL99j , 
used in [DeL03, Dc L04[ lDeL06j to generate approximations of surfaces' gener- 
ality stability zones against their analytical boundaries. 

Indeed no simple algorithm to generate the analytical equations of the bound- 
aries of stability zones for a generic function is known so far but we know that all 
directions belonging to the same stability zone share the same soul and therefore 
it is possible to get an approximate picture of the set of generalized stability 
zones by evaluating the soul of some (big) set of rational directions. For ex- 
ample in all cases examined so far, thanks to the high level of symmetry, we 
could limit our analysis to the directions contained in the the triangle with sides 
(0 : : 1), (1 : 1 : 1) and (1:0:1) (concretely, to all directions {m : n : N), 
< n < m < N, for N = 10 2 , 10 3 , 10 4 ). 
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Figure 9: (a) Picture of the fractal in the square [0, l] 2 in the z — 1 chart of RP 2 
obtained by evaluating the soul of every stability zone in the lattice (n, m, N), n,m = 
1,...,10 3 , N = 10 3 ; (b) the same picture compared with the analytical boundaries of 
the 797161 stability zones of the first 12 levels of recursion. Out of the 961367 rational 
directions for which a label was numerically found by the algorithm, 455654 belong 
to those 797161 stability zones and all of them turn out to lie within the analytical 
boundaries of the corresponding zone. 

Note that rational directions in this setting are of paramount importance 
because their (non-critical) leaves are compact (in T 3 ) and therefore can be in 
principle approximated with error as small as wished and therefore the corre- 
sponding soul can be, at will, evaluated exactly through numerical calculations. 
Moreover, rational directions are dense in every stability zone and therefore (in 
principle) we do not lose any picture detail by restricting our analysis to them. 

Below we present the algorithm we use to retrieve the soul (if any) of a 
rational direction (m : n : N) in this particular case, where we have all saddles 
of "monkey" type: 

NO choose a representative &j, i = 1,2,3, for the cycles on C which are respec- 
tively homologous in T 3 to the x, y and z axes; 

Nl retrieve the intersection between C and the plane 

m(x - 1/4) + n(y - 1/4) + N(z - 1/4) = ; 

N2 follow the three critical loops and, if no saddle connection is detected, 
store them in the variables 01,2,3, otherwise exit; 

N3 evaluate the homology class of 01,2,3 m T 3 ; 

N4 if exactly one of the three loops is non-zero in T 3 then evaluate its inter- 
section numbers with the {bi} and set this triple as the soul corresponding 
to the direction (m,n,N), otherwise exit. 
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The result of sampling the triangle T with a 10 3 x 10 3 lattice are shown in 
Figure[9]and turn out to be in perfect agreement with the analytical boundaries. 
An evaluation of the fractal box dimension based on these numerical data leads 
to a value of about 1.7, which is also in very good agreement with the evaluation 
obtained from of the analytical boundaries made in sec. 13.21 
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